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Abstract- This  paper  presents  a  nonlinear  extension  of  the  affine 
transform  as  used  in  the  iterated  function  systems  (IFS)  to  com¬ 
press  signals  and  data.  Compared  with  0ien  and  Narstad’s 
orthogonal  transform  with  compression  ratio  of  6.0:1  for  the 
electrocardiogram  (ECG)  signal,  the  nonlinear  approach  pro¬ 
duces  higher  compression  ratios  of  10.2:1,  and  is  more  flexible 
in  finding  the  corresponding  strange  attractor  of  the  ECG  sig¬ 
nal.  It  can  model  the  QRS  complex  of  the  ECG  signal  well, 
which  has  been  a  problem  for  the  affine  transform  in  fractal 
compression. 

Keywords-ECG,  strange  attractor,  IFS,  affine  transform. 

I.  Introduction 

Data  compression  based  on  fractals  commenced  about 
ten  years  ago  [1],  More  specifically,  the  iterated  function  sys¬ 
tems  (IFS)  approach  to  data  compression  has  been  developed 
based  on  the  fundamental  property  of  fractal  objects:  self¬ 
similarity.  Thus,  our  discussion  begins  from  a  fractal  signal. 
A  random  discrete  process  x(n)  defined  for  all  integers, 
-°o  <  n  <  °° ,  is  said  to  be  statistically  self-similar  if  its  statis¬ 
tics  are  invariant  to  its  dilations  and  contractions.  Thus,  x(n) 
is  statistically  self-similar  with  a  parameter  D  if  for  any  real 
a  >  0 ,  it  obeys  the  following  scaling  relation  [2] 

P  -D 

x(n)  =  a  x(a  n)  (1) 

P 

where  =  denotes  equality  in  a  statistical  sense. 

To  employ  the  multiplicative  (rather  than  additive) 
invariance  of  the  dilation  and  contraction  of  the  fractal  sig¬ 
nal,  Barnsley  proposed  a  fractal  compression  technique,  the 
IFS,  which  uses  an  affine  transform  to  map  the  process  onto 
itself  [1],  Such  an  affine  transform,  W,  for  a  one-dimensional 
process  x(n)  is  defined  as 
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where  a,  c,  d,  e,  and/are  the  transform  coefficients. 

The  IFS  compression  technique  partitions  the  process 


into  domain  blocks  and  range  blocks,  and  finds  a  set  of  affine 
transforms  between  the  domain  blocks  and  range  blocks. 
Computational  complexity,  compression  ratio,  reconstruc¬ 
tion  quality  of  the  signal,  and  convergence  are  important 
issues  in  the  IFS  approach  [1],  All  these  issues  are  related  to 
the  quality  of  the  matching  between  the  blocks  which  results 
in  the  affine  transform.  The  role  of  the  affine  transform  is  to 
fit  the  range  blocks  with  the  domain  blocks  optimally.  The 
compression  ratio  and  reconstruction  quality  are  decided  by 
the  optimization  algorithm  of  the  affine  transform.  The  cur¬ 
rent  transform  in  the  IFS  described  by  (2)  is  a  linear  affine 
transform.  Its  fitting  capability  is  not  satisfactory  for  signals 
such  as  image,  speech,  and  electrocardiogram  (ECG),  which 
have  nonlinear  characteristics.  In  addition,  the  original  IFS 
search  algorithm  is  very  time  consuming  due  to  the  direct 
comparison  of  a  large  number  of  domain  blocks  [3],  An 
order  of  magnitude  improvement  in  the  search  time  was 
achieved  through  a  frequency-sensitive  neural-network 
search  [4]  [5].  It  is  necessary  to  find  another  model  for  the 
IFS  to  improve  compression  performance  and  reduce  the 
search  time. 

A  more  flexible  way  to  model  complicated  signals  is  to 
use  a  nonlinear  transform.  We  use  a  polynomial  expansion  to 
model  the  ECG  signal  because  it  is  an  efficient  way  to  model 
arbitrary  functions.  Convergence  is  an  important  issue  in  the 
signal  reconstruction  of  the  IFS  for  the  proposed  transform. 
Although  we  have  not  proved  it  in  theory,  our  experiments 
show  that  the  new  extended  affine  transform  gives  conver¬ 
gent  results.  This  nonlinear  model  also  demonstrates  the 
flexibility  in  modelling  various  complexities  in  the  ECG  sig¬ 
nal.  We  shall  show  that  it  can  give  much  better  reconstruc¬ 
tion  quality  under  the  same  compression  ratio,  as  compared 
with  a  corresponding  linear  model.  The  difficulty  in  model¬ 
ling  the  QRS  complex  in  ECG  signals  [6]  is  solved  with  the 
proposed  technique.  The  new  transform  combined  with  com¬ 
plexity  measure  can  reduce  the  computational  complexity 
[7].  Experimental  results  are  presented  and  compared  with 
the  standard  IFS  and  other  techniques. 

For  clarity  and  completeness,  a  continuous  ECG  signal 
is  first  obtained  from  multi-channel  sensors,  then  amplified 
and  filtered  to  satisfy  the  Nyquist  sampling  theorem,  then 


Report  Documentation  Page 

Report  Date  Report  Type 

25  Oct  2001  N/A 

Dates  Covered  (from..,  to) 

Title  and  Subtitle 

Contract  Number 

A  New  Nonlinear  Transform  for  IFS  Compression  of  ECG  and 
Other  Signals 

Grant  Number 

Program  Element  Number 

Author(s) 

Project  Number 

Task  Number 


Work  Unit  Number 

Performing  Organization  Name(s)  and  Address(es)  Performing  Organization  Report  Number 

Department  of  Electrical  and  Computer  Engineering  Signal  and 
Data  Compression  Laboratory  University  of  Manitoba, 

Winnipeg,  Manitoba,  R3T  5V6,  Canada 

Sponsoring/Monitoring  Agency  Name(s)  and  Address(es)  Sponsor/Monitor’s  Acronym(s) 

US  Army  Research,  Development  &  Standardization  Group  - 

(UK)  PSC  802  Box  15  FPO  AE  09499-1500  Sponsor/Monitor’s  Report  Number(s) 

Distribution/ Availability  Statement 

Approved  for  public  release,  distribution  unlimited 

Supplementary  Notes 

Papers  from  23rd  Annual  International  Conference  of  the  IEEE  Engineering  in  Medicine  and  Biology  Society,  October 
25-28,  2001,  held  in  Istanbul,  Turkey.  See  also  ADM001351  for  entire  conference  on  cd-rom..  The  original  document 
contains  color  images. 

Abstract 


Subject  Terms 


Number  of  Pages 

5 


2  of  5 


sampled  to  form  a  discrete  signal,  and  finally  quantized  to  strange  attractor  of  the  object,  it  is  reasonable  to  assume  [1] 
form  a  digital  ECG  signal.  This  paper  is  concerned  with  the 

digital  ECG  form  only.  g(x(n {),  n {)  =  x{n2)  (5) 

n.  The  Nonlinear  IFS  (NIFS)  where  x(n  | )  is  mapped  to  x(n0)  under  the  condition 


The  IFS  data  compression  is  realized  through  finding  a 
set  of  affine  transforms  between  range  pool  and  domain 
pool.  The  role  of  the  affine  transform  is  to  fit  specific  range 
blocks  with  the  domain  blocks  optimally.  The  compression 
performance  is  directly  decided  by  the  capability  of  the  aff¬ 
ine  transform  that  models  the  signal.  The  current  transform 
used  in  the  IFS  is  a  linear  affine  transform.  It  may  model  per¬ 
fectly  some  signals  which  have  strict  inherent  linear  relation¬ 
ship  (such  as  the  Cantor  set,  Koch  curve,  Sierpinski  carpet, 
and  Julia  set)  [8],  Unfortunately,  many  signals  in  real  life 
have  relationships  that  are  more  complicated  than  linear.  To 
model  such  signals,  the  linear  affine  transform  is  inadequate. 
Image  compression  by  the  IFS  fractal  method  is  an  example 
of  how  difficult  it  is  to  find  a  suitable  mapping  between  a 
range  block  and  the  domain  blocks.  A  complicated  scheme 
that  requires  a  point-by-point  search  combined  with  a  varia¬ 
ble  size  of  the  domain  blocks  is  necessary  to  find  a  reasona¬ 
ble  mapping,  even  though  it  may  fail  occasionally.  This 
transform  causes  two  problems:  (i)  poor  compression  per¬ 
formance  and  (ii)  a  very  time-consuming  search. 

Instead,  a  generalized  transform  which  may  represent  an 
arbitrary  function  is  needed  to  model  such  complicated  sig¬ 
nals.  For  a  one-dimensional  signal  x(n),  a  generalized  trans¬ 
form  is  defined  as 


w( n']  =  an  +  e 

w  g(x,  n)_ 


(3) 


an^  +  e  =  n2  (6) 

Equations  (4),  (5),  and  (6)  constitute  our  nonlinear  con¬ 
traction  transform  for  finding  the  strange  attractor  of  the 
object.  When  the  expansion  order  k  is  equal  to  1,  the  nonlin¬ 
ear  transform  becomes  the  traditional  affine  transform. 
Therefore,  one  sees  that  the  nonlinear  transform  may  lead  to 
more  flexible  mappings  between  the  range  blocks  and  the 
domain  blocks  than  the  linear  affine  transform.  Such  flexibil¬ 
ity  may  result  in  two  advantages:  (i)  improving  compression 
performance  and  (ii)  speeding  up  the  search. 

We  shall  now  discuss  how  to  use  the  least  square  error 
technique  to  determine  the  coefficients  in  (4)  and  (6).  For  a 
lossless  fractal  compression  technique,  a  transform  must  be 
found  to  satisfy  (5)  with  a  measured  signal  exactly.  In  prac¬ 
tice,  it  is  reasonable  to  use  an  approximation  to  replace  the 
strict  equality  condition  between  g(x,  n)  and  x{n)  in  the 
lossy  fractal  compression.  Then  there  is  an  error 

ri 

E-  Y  [gWnj),  «j) -x(«2)]“  (7) 

n2  =  1 

where  r,  is  the  length  of  the  range  block. 

A  set  of  optimal  transform  coefficients  can  be  obtained 
by  applying  the  least  square  error  approach  on  (7). 

III.  Experimental  Results  and  Analysis 


where  the  coefficient  a  is  limited  to  the  range  [0,  1 )  to  guar¬ 
antee  contraction  of  signals  in  time.  If  g(x,  n)  is  also  a  con¬ 
tractive  transform,  then  (3)  satisfies  the  convergence 
condition  of  the  IFS. 

To  find  an  appropriate  g(x,  n) ,  it  is  not  sufficient  to  use 
the  convergence  condition  only.  According  to  Taylor’s 
remainder  formula,  a  function  may  be  expanded  as  a  kth 
order  polynomial  format  if  its  first  k  -  1  derivatives  are  con¬ 
tinuous  and  its  kth  derivative  exists  [9],  Thus,  a  kth-order 
polynomial  is  used  to  approximate  the  g{x,n) 

k  p 

g(x,  n)  ~  Y  Vk,  p,qeZ  (4) 

p  =  0q  =  0 

where  coefficients  b  is  decided  by  the  function  g(x,  n) 

pq 

according  to  Taylor’s  remainder  formula.  Since  g(x,n)  is 
unknown,  it  may  be  determined  with  the  inverse  procedure 
by  finding  the  b  .  The  definition  of  the  strange  attractor  of 
objects  gives  a  way  of  how  to  find  the  b .  To  obtain  the 


The  new  nonlinear  IFS  transform  is  applied  to  compress 
a  one-dimensional  ECG  signal.  To  compare  compression 
performance,  the  traditional  IFS  technique  is  also  imple¬ 
mented.  The  performance  of  a  compression  algorithm  is 
dependent  mainly  on  two  parameters:  compression  ratio  and 
reconstruction  (distortion)  error.  The  compression  ratio,  Rcr, 
is  defined  as  [10] 

N 

Y  size^XorM^ 

Rcr  =  ^ -  (8) 

Y  size^xcmM^ 

n  =  1 

where  x  An)  and  x  .(/;)  represent  the  original  and  the 
compressed  signal,  respectively.  Their  corresponding  sample 
lengths  are  N  and  N  l . 

A  common  distortion  measure  is  the  following  percent 
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root-mean-square  difference  (PRD)  to  measure  recon¬ 
structed  error  of  signals  [11] 


PRD 


N 

E 

[xoW(n)-xrec(n)l2 

n  =  1 

N 

Y  Xnrhn ) 

X100%  (9) 


n  =  1 


where  xr  (n)  represents  the  reconstructed  signal.  In  the 
processing  of  reconstruction  error,  it  is  found  that  the  PRD 
will  change  when  the  signals  shift  along  amplitude  direction. 
A  normalized  PRD  measure  is  proposed  in  [12]  for  ECG  sig¬ 
nals 


NPRD  = 


i  y  r  xori(n)-xreM)  l2 
AN  Lf  max(xori)  -  min(xori ) 

V  n  =  1 


x  100%  go) 


Notice  that  although  (9)  is  not  a  good  error  measure  in 
ECG  data  compression,  it  is  used  in  this  paper  to  measure 
reconstruction  error  for  convenience  to  compare  our  results 
with  other  research. 


Coefficient  Amplitude 

Fig.  1  Distribution  of  coefficients  (a)  Z?10  and  (b)  b2 0  ■ 

Prior  to  the  calculation  of  the  compression  ratio  and 
reconstruction  error,  the  coefficients  in  g(x,  n)  must  first  be 
quantized.  Figure  1  shows  two  distributions  of  coefficients 


b  |Q  and  b1{]  with  the  ECG  signal  sampled  at  360  Hz  with 
1 1  bit  resolution,  as  taken  from  the  MIT-BIH  ECG  database 
[13],  We  do  experiment  using  the  ECG  data  contained  in  the 
file,  x_100.txt,  which  has  a  10-minute  recording  of  the  ECG 
signal.  The  Lloyd-Max  nonuniform  quantizer  is  used  to 
quantize  the  coefficients  because  their  distribution  is  close  to 
a  Laplacian  [14]. 

To  improve  compression  performance,  the  coefficient  a 

—171 

is  set  to  2  where  m  =  1,  2,  3  ,  to  increase  the  chance  for 
finding  the  optimal  mapping.  The  length  of  range  blocks  is 
also  a  variable,  changing  according  to  2  with  the  integer  / 
varying  from  2  to  6.  Figures  2(a)  and  (b)  give  the  size  distri¬ 
butions  of  range  blocks  for  the  traditional  and  the  nonlinear 
IFS  (NIFS)  approach,  respectively.  The  NIFS  can  use  the 
most  frequent  range  block  size  of  32  sample  points  to  model 
the  ECG  signal  compared  with  4  sample  points  used  by  the 
linear  affine  transform.  The  size  distributions  of  range  blocks 
show  that  more  large  range  blocks  are  found  by  the  NIFS 
than  by  the  traditional  IFS,  thus  leading  to  an  improved  com¬ 
pression  ratio.  In  our  experiment,  the  mean  range  size  is 
about  30.0  for  the  NIFS  and  26.9  for  the  IFS.  It  demonstrates 
that  the  nonlinear  transform  can  model  signals  more  flexibly. 


Range  Block  Size 


Fig.  2  Range  block  size  distributions  for  (a)  the  traditional  affine  transform 
and  (b)  the  nonlinear  transform. 


Figures  3(a)  and  (b)  give  address  distribution  of  the  opti¬ 
mally  mapped  domain  blocks  found  by  the  traditional  IFS 
and  the  NIFS  approach,  respectively.  The  address  distribu¬ 
tion  of  the  mapped  blocks  given  by  the  NIFS  demonstrates 
more  concentrated  distribution  around  the  first  point  of  the 
search.  Compared  with  full  distribution  in  the  segment  of 
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1024  points  of  the  IFS,  the  NIFS  gives  distribution  of 
domain  block  address  within  110  point  range.  It  shows  that 
the  optimal  mapping  can  found  by  searching  fewer  domain 
blocks.  This  can  be  useful  in  developing  a  fast  search  algo¬ 
rithm  in  the  future. 


Domain  Address 


Fig.  3  Address  distribution  of  the  optimally  mapped  domain  blocks  found 
by  (a)  the  traditional  and  (b)  the  NIFS  approach. 

Figures  4(a),  (b),  and  (c)  show  the  original  ECG  signal, 
its  reconstruction,  and  the  reconstructed  error.  From  these 
figures,  one  can  see  that  the  NIFS  compresses  and  recon¬ 
structs  the  ECG  signal  extremely  well.  The  new  technique 
has  the  remarkable  property  of  modelling  the  QRS  complex 
of  the  ECG  signal  almost  perfectly,  which  is  a  difficult  issue 
for  the  linear  affine  transform  and  the  orthogonal  transform 
[6].  The  reconstruction  error  is  limited  to  a  small  range  only. 
It  is  mainly  due  to  the  noise  in  the  signal. 

Table  1  shows  the  influence  of  the  order  k  of  the  Taylor 
series  expansion  and  quantization  resolution  on  the  compres¬ 
sion  performance  of  the  ECG  signal.  By  changing  the  expan¬ 
sion  order  from  1st  to  4th,  as  well  as  the  resolution  of  the 
quantizer  from  7  bits  to  14  bits,  various  PRD  and  Rcr  are 
obtained.  By  monitoring  the  changes  among  them,  optimal 
parameters  for  the  NIFS  on  the  ECG  signal  can  be  found. 

Table  1  shows  that  when  the  resolution  increases,  the 
Rcr  increases.  The  PRD  is  controlled  by  the  error  threshold. 
When  quantizer  takes  more  than  11  bits,  the  PRD  and  the 
Rcr  change  little  with  the  resolution  increasing  further.  In  the 
future,  we  will  use  different  quantizer  for  different  coeffi¬ 
cient. 


0.5  1  1.5  2  2.5 

Reconstructed  Error  [s] 


Fig.  4  The  ECG  signal  compression:  (a)  original  ECG  signal, 
(b)  reconstructed  signal,  and  (c)  reconstructed  error. 


Table  1 :  Compression  performance  investigation  of  the  NIFS  with 
quantization  resolution  and  polynomial  expansion  change  on  the  ECG  signal. 


Resolution 

(bit) 

k=l 

k=2 

k=3 

k=4 

PRD 

Rcr 

PRD 

Rcr 

PRD 

Rcr 

PRD 

Rcr 

7 

6.0% 

9.1 

5.8% 

6.0 

5.6% 

3.0 

5.5% 

2.2 

8 

5.8% 

9.0 

5.4% 

5.4 

5.7% 

4.9 

5.7% 

4.0 

9 

5.8% 

9.5 

5.7% 

8.0 

5.6% 

6.1 

5.8% 

5.3 

10 

5.7% 

10.5 

5.6% 

7.6 

5.3% 

6.0 

5.8% 

6.1 

11 

5.8% 

10.6 

5.9% 

9.9 

5.8% 

7.9 

5.9% 

7.4 

12 

5.8% 

10.7 

5.8% 

9.9 

5.9% 

8.6 

5.8% 

7.7 

13 

5.9% 

10.6 

5.8% 

10.2 

6.0% 

9.2 

6.0% 

8.3 

14 

5.8% 

11.1 

5.8% 

10.2 

5.7% 

9.2 

5.4% 

8.4 

Threshold 

0.0115 

0.0117 

0.0123 

0.0130 

One  may  predict  that  the  expansion  order  increase  will 
result  in  a  PRD  decrease  because  of  the  more  powerful  mod¬ 
elling  ability  of  the  NIFS.  However,  a  higher  order  results  in 
more  coefficients,  thus  leading  to  a  lower  compression  ratio. 
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Consequently,  there  must  be  a  trade-off  between  the  RMS 
and  the  Rcr  for  choosing  the  optimal  order  k.  Although  Table 
1  demonstrates  that  the  change  of  expansion  order  has 
almost  no  influence  on  the  PRD  (because  the  PRD  is  control¬ 
led  by  the  error  threshold),  it  results  in  different  Rcr.  The 
expansion  order  of  2  may  compete  with  the  compression 
ratio  achieved  by  the  linear  IFS. 

It  was  found  that  the  NIFS  achieves  almost  the  same 
compression  ratio  as  the  traditional  IFS  under  the  same 
reconstruction  error  in  the  ECG  signal  compression.  When 
PRD  =  5.8%,  0ien  and  Narstad  achieved  6.0:1  with  their 
IFS  ECG  compression  by  orthogonal  transform  [6].  With  the 
optimal  choice  of  order  2  and  14  bit  nonuniform  quantizer, 
we  get  the  compression  ratio  about  11.1:1  for  the  affine 
transform  and  about  10.2:1  for  the  NIFS. 

IV.  Conclusion 

In  this  paper,  we  use  a  nonlinear  IFS  (NIFS)  transform 
to  compress  a  one-dimensional  ECG  signal.  The  first  impor¬ 
tant  issue  in  this  technique  is  the  convergence  of  the  strange 
attractor  reconstruction.  We  have  seen  that  the  signal  recon¬ 
struction  is  successful.  The  experiment  shows  that  the  non¬ 
linear  approach  may  compete  with  the  traditional  IFS 
technique  for  the  ECG  signal  compression.  With  the  optimal 
choice  of  order  and  nonuniform  quantizer,  the  extended 
transform  achieves  a  maximal  compression  ratio  of  10.2:1, 
which  is  higher  than  that  of  6.0:1  obtained  by  0ien  and 
Narstad  under  the  same  PRD. 

Another  advantage  of  the  new  extended  transform  is  that 
it  can  model  the  QRS  complex  of  the  ECG  signal  very  well, 
which  has  been  a  problem  for  the  affine  transform  and  the 
orthogonal  transform  in  fractal  compression.  Fewer  blocks 
are  used  to  model  the  ECG  signal  by  the  NIFS  than  that  by 
the  linear  IFS.  It  reveals  that  the  NIFS  technique  has  a  more 
powerful  modelling  ability. 

The  address  distribution  of  the  optimally  mapped 
domain  blocks  around  the  first  point  in  the  search  in  Fig.  3(b) 
has  important  significance.  It  can  be  employed  to  speed  up 
the  search  of  finding  optimal  mappings  between  the  range 
blocks  and  domain  blocks.  In  another  paper,  we  describe 
how  to  use  the  extended  transform  to  develop  a  fast  NIFS 
algorithm. 
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